Three body resonances in two meson-one baryon systems 
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We report four E's and three A's, in the 1500 - 1800 MeV region, as two meson - one baryon S-wave 
(l/2) + resonances. We solve Faddeev equations in the coupled channel approach. The invariant 
mass of one of the meson-baryon pairs and that of the three particles have been varied and peaks 
in the squared three body T- matrix have been found very close to the existing S = -1, J p = l/2 + 
low lying baryon resonances. The input two-body t-matrices for meson-meson and meson-baryon 
interaction have been calculated by solving the Bethe-Salpeter equation with the potentials obtained 
in the chiral unitary approach. 

PACS numbers: 21.45.+V, 25.75.Dw 



Hyperon physics is still relatively unexplored and there 
are many open problems to be studied. The hyperon res- 
onances, for example, have a poor status as compared 
to the nucleon ones |l[ . Though some of them, like the 
low lying J p = 1/2" resonances, A(1405), A(1670) 
can be represented as dynamically generated and rela- 
tively well understood states within the unitary chiral 
models, the low lying J p = l/2 + domain remains far 
less understood, both experimentally and theoretically. 
For instance, quark models seem to face difficulties in 
reproducing properties of the resonances in this sector 
0. The neat reproduction of the low lying 1/2 - states 
in the S-wave meson-baryon interaction, using chiral dy- 
namics, suggests that the addition of a pseudoscalar me- 
son in S-wave could lead to an important component of 
the structure of the l/2 + resonances. There exist re- 
sults which hint towards this possibility, e.g., in Q it 
was found that the two meson cloud gave a sizable con- 
tribution to the mass in the spectrum of the l/2 + baryon 
antidecuplet. Chiral dynamics has been used earlier in 
the context of the three nucleon problems, e.g., in [J]. 
In this article we present the first study of two meson - 
one baryon systems where chiral dynamics is applied to 
solve the Faddeev equations. As shall be described be- 
low, our calculations for the total S = -1 reveal peaks 
in the wKN system and its coupled channels which we 
identify with the resonances £(1770), £(1660), £(1620), 
£(1560), A(1810) and A(1600). 

We start by taking all combinations of a pseudoscalar 
meson of the 0~ SU(3) octet and a baryon of the l/2 + 
octet which couple to S = — 1 with any charge. For some 
quantum numbers, the interaction of this two body sys- 
tem is strongly attractive and responsible for the gen- 
eration of the two A(1405) states @ and other S = 
-1 resonances. We shall assume that this two body 
system formed by KN and coupled channels remains 
highly correlated when a third particle is added, in the 
present case a pion. Altogether, we get twenty- two 



coupled channels for the net charge zero configuration: 



7r°7r £°, 



-£+ 



7r°77£°, tt^A, tt°K+E-, 7r°iv°S , n+K-n, 7r+7r°£-, 
7r+7r-£°, tt+tt-A, 7T+77£-, ir+K°E-, n-K°p, 7r-7r°£+, 
7r~7T + £°, 7r~7r + A, 7r~?7£ + , 7r~ K + E° . We assume the 
correlated pair to have a certain invariant mass, y / S23, 
and the three body T-matrix is evaluated as a function 
of this mass and the total energy of the three body sys- 
tem. 

The input required to solve the Faddeev equations, i.e., 
the two body t-matrices for the meson-meson and meson- 
baryon interactions have been calculated by taking the 
lowest order chiral Lagrangian following @, 0, H, 0] and 
using the dimensional regularization of the loops as done 
in [1, [l(| , where a good reproduction of scattering am- 
plitudes and resonance properties was found. Instead, 
a cut off could also be used to regularize the loops as 
shown in 0, H3] • Improvements introducing higher order 
Lagrangians have been done recently, including a theo- 
retical error analysis [ll[ which allows one to see that 
the results with the lowest order Lagrangian fit perfectly 
within the theoretical allowed bands. 

A shared feature of the recent unitary chiral dynam- 
ical calculations is the on-shell factorization of the po- 
tential and the t-matrix in the Bethe-Salpeter equation 
d 0, M, HI HUH, which is justified by the use of 
the N/D method and dispersion relations [13, EH- Al- 
ternatively, one can see that the off-shell contributions 
can be reabsorbed into renormalization of the lower or- 
der terms @, 0| ■ We develop here a similar approach for 
the Faddeev equations. 

The full three-body T-matrix can be written as a sum 
of the auxiliary T-matrices T 1 , T 2 and T 3 [l7| 



(1) 



where T l , i = 1, 2, 3, are the normal Faddeev partitions, 
which include all the possible interactions contributing 
to the three-body T-matrix with the particle i being a 
spectator in the last interaction. The Faddeev partitions 
satisfy the equations 
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T' 



t i S 3 (k{ - h) + fg ij T j + t l g ik T k 



(2) 
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where ki (fc/) is the initial (final) momentum of the ith 
particle in the global center of mass system, t % is the two- 
body t-matrix for the interaction of the pair (jk) and g % i 
is the three-body propagator or Green's function, with 
j + k ? i = 1, 2, 3 

Iterating Eq.([2]) and removing the term with 5 3 (k- — 
fcj), which correspond to a "disconnected diagram" [H| . 
will give 

v = /'//'•'/•• + e g ik t k + t i g ij t j g jk t k + iy^v'v + 

+fg ik t k g k H j + e g lk t k g kt e + ■■■ (3) 

The first two terms of the Eq.Q, for the case i=l, 
are represented diagrammatically in Fig[I] where the t- 
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FIG. 1: Diagrammatic representation of the terms (a) t^g 12 ^ 

(b) ty i3 >t 3 . 

matrices are required to be off-shell. However, the chiral 
amplitudes, which we use, can be split into an "on-shell" 
part (obtained when the only propagating particle of the 
diagrams, labeled with £?, nt in FigfTJ is placed on-shell), 
which depends only on the cm energy of the interacting 
pair, and an off-shell part proportional to the inverse of 
the propagator of the off-shell particle. This term would 
cancel the particle propagator (for example that of the 
3rd particle in the Figfjji) resulting into a three body 
force (Figj2^i) . In addition to this, three body forces also 
stem directly from the chiral Lagrangians (FigOj). 





FIG. 2: The origin of the three body forces (a) due to cancel- 
lation of the propagator in Fig[]Ja) with the off-shell part of 
the chiral amplitude, (b) at the tree level. 

We find that the sum of the off-shell parts of all the 
six t % g l Hi terms, together with the contribution from 
FiglUb) cancels exactly if the SU(3) limit is considered 
and the momentum transfer for the baryon is assumed 
to be small. In a realistic case we find this sum to be 
smaller than 5% of the total on-shell contribution. Hence, 
only the on-shell part of the two body (chiral) i-matrices 
is significant. The diagrams in Fig.l can then be ex- 
pressed mathematically (reading the diagrams from right 
to left as a convention) as i 1 ^ 12 ^ 2 and t 1 <7 13 t 3 , respec- 
tively, where the ^'-matrices depend only on the center 
of mass energy of the interacting particles. 



The t l g % Hi terms correspond to the situation where 
there are no loops and hence the g* J propagators are writ- 
ten in terms of the on-shell variables 



1 




Eiik'i) - Ejikj) - E k (k[ +kj) + ie 



with t/s being the total energy in the global CM system. 



Ei = y kf + mf is the total energy of the particle I and 
Ni is a normalization constant (TV; = 1 for the meson- 
meson interaction and Ni = 2Mi for the meson-baryon 
interaction, where M; is the corresponding baryon mass) 
and D is the number of particles propagating between 
two consecutive interactions. 

These propagators (and all other angle dependent ex- 
pressions in the formalism) are projected in S-wave. A 
proper Lorentz boost has been made for transformation 
of the momenta from the center of mass frame of two par- 
ticles to the global center of mass frame whenever needed. 
A technical remark is here in order: to avoid the evident 
divergence in these on-shell propagators at the threshold 
of a channel and to continue the Faddeev equations an- 
alytically below the threshold, we fix the momentum of 
the particle to a minimum value, P m i n = (for example, 
50 MeV). It should be mentioned that the results are al- 
most insensitive to a change in the value of the P m in, 
since a change in P m i n of ~ 40- 50% results into a shift 
in the peak position by less than 2 MeV. 

The first term with a non trivial structure, from the 
point of view of the on-shell factorization of the t- 
matrices in the Faddeev equations, is the one involving 
three successive pair interactions, where a loop function 
of three particle propagators appears for the first time. 
We show the diagrams with such a structure for the T 1 
partition in Fig[3]Ja-d). 
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FIG. 3: Different diagrams involving three pair interactions 
corresponding to the T 1 partition. 

We write all these t l g l H^g^ k t k terms as 
i*G*' fc (iV fc t fc )U-./.di where 



Qij k 



d 3 k" Ni N„ 



F ijk (^/s,k") 



(2tt)3 2E l 2E m ^ - El (k") - E m (k") + ie 

(4) 
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with i 7^ j , j =/= k, i =/= I =/= 77i and yjsi rn is the invariant 
mass of the (lm) pair. Eq.Q consists of the two particle 
propagator in the first cut (shown as a dashed line in 
Fig. 3 (a) as an example) and the F l3k function, which is 
defined as 



F ljk = t 3 (y 



9 jk \off- 



shell 



9 jk \on- 



shell 



(5) 

Sint{k") = s — m'j — 2y/sEj(k") denotes the invariant 
mass required to calculate the t 3 -matrix in the loop in- 
tegral of G i3 K The term g 3k \ m l shell [t 3 (^(kr))}- 1 
appearing in F ljk of Eq.([S]) cancels the {P g 3k )\ on -shell 
of VG^figH^onsheu and produces the t i g ij t 3 g 3k t k 
term with g l H 3 depending on the proper off-shell vari- 
able, k" of the loop. This procedure allows us to render 
the Faddeev equations into a set of algebraic equations, 
as we see below. To regularize the integrals of Eq. we 
shall use the cut off as in [(| 0] which is of the order of 
1 GeV. The results are rather insensitive to this cut off, 
due to the convergence of these loops which involve three 
propagators. 

So far we have discussed diagrams with two or three 
t-matrices. It has been shown that the introduction of 
a third interaction replaces the propagator g by a loop 
function G. The formalism is further developed by mak- 
ing the same substitution whenever a new interaction is 
added. 

The Faddeev partitions of Eq.([3]) in this prescription 
can be re-written as 



/ = \r<ri* ■ r(r^ r ■ re""/-,,"/ 

+ {fg lk t k +t l G tk n k g k H 1 



} 



t i G iki t k g ki f + 



— 1 R 



T, 



"ik 



(6) 



where the T l partition has been rewritten in terms of the 
two infinite series T% and T lk , which sum all the dia- 
grams with the last two interactions written in terms of 
t 1 , P and t l , t k , respectively. Hence, we obtain six par- 
titions (which is double of those in the original Faddeev 
equations); 



H 3 + f 



G ijl T 3 ^ 



Qijkrpjk 



(7) 



with i j =/= k. 

It remains to define the invariant masses on which the 
two body t-matrices and the propagators depend. The 
expression for the s±2 and S13 obtained from the energy 
conservation in terms of the external (on-shell) variables 
is 



Si, = S + 771 1. — 



2 ~ Ex) (s 2 3 + ml - m) ) 



S23 



(8) 



with E\ — (s — S23 + mf)/(2y / i). However, it should 
be noted that the t-matrices t 3 {y/Si n t{k")) in the loop 
(Eq.Q) are calculated in terms of a running variable as 
required. 



There are two independent variables in the formalism 
y/s, \/s23, as a function of which we plot the squared 
Tj^-matrix (T£ = £)(Tft ~ ) ), since the t l g l3 t 3 

terms evidently do not have a resonance structure and 
just provide a background to the amplitudes. 

|T* R ] 2 (1CT 9 MeV 6 ) 



< s 23 (MeV) 




15go 1620 1650 1680 



1470 1500 1530 1560 

4 s (MeV) 



FIG. 4: Two E resonances in the 7r7rE amplitude in I = 1, 1^ 
= 2 configuration. 

We now report the four isospin one states found in our 
study. In Figdl we show a plot of the squared T^-matrix 
and its projection, for 7T7tE — > 7T7tE in the total isospin 
1=1 configuration obtained by keeping the two pions 
in isospin I„ = 2. We see two peaks; one at yfs = 1656 
MeV with the full width at half maximum ~ 30 MeV and 
another at y/s = 1630 MeV with T = 39 MeV. We identify 
the peak at y/s = 1656 MeV with the well established 
E(1660 — il00/2) [l[ as a resonance in the 7T7tE system, 
which is a new finding. It is interesting to recall that 
the excitation of this resonance is claimed in the study 
of the K~ p -> 7r 7T°E reaction (20|. Note that since 
we are plotting the squared amplitude, which would be 
proportional to a cross section of a certain process, we 
can associate our results to the ordinary masses of the 
PDG and not to the "pole positions" also quoted there. 

The peak in squared T^-matrix observed at 1630 MeV 
with a width of 39 MeV needs a special attention. The 
two-star resonance £(1620) Q, though listed as a l/2~ 
state, seems to be a very unclear case. The partial wave 
analysis and the production experiments have been kept 
separately in [l[ since it is difficult to know the quantum 
numbers from the production experiments and if more 
than one resonance contributes to a single bump. Inter- 
estingly, there is a l/2 + state found by the partial wave 
analysis work of Martin et. al. [2l[ in this region. The 
authors of i2l"j use a multichannel partial wave analysis 
of the KN data and find a resonance at 1597 MeV. This 
result has however been listed under the £(1660) in [l|. 
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FIG. 5: A A resonance in the -kKN amplitude at 1568 MeV 
in / = 0, I nR = 1/2. 

Another partial wave analysis of the KN — * An reac- 
tion made by Armenteros et. al. [13] find a 1/2+ Pn 
resonance at 1610 MeV with a width of 60 MeV. These 
findings would provide some phenomenological support 
to our claim of a l/2 + £ resonance around ~ 1620 MeV. 

We find two more peaks in the I — 1 sector; one at 
sfs = 1590 MeV with a width ~ 70 MeV in I = 1, 1 % 
= state and another at yfs = 1790 MeV with T = 24 
MeV in 7 = 1, Z,r = 2 case. The former one supports the 
existence of the £(1560) "bump", whose spin-parity is 
unknown [l[. Our results would associate a l/2 + to the 
spin-parity of this resonance. The latter finding supports 
the one-star £(1770). 

Next, we discuss the three isospin zero states obtained 
in these calculations. First we look at states observed in 
■kKN with I^g = 1/2. Two peaks in the A(1600) MeV 
region have been found at yfs — 1568 with a width of 60 
MeV (which is shown in FigO and at 1700 MeV with 
T=136 MeV. One should note that the PDG quotes a 
mass for the A(1600) between 1560 MeV and 1700 MeV 
and the width between 50 MeV and 250 MeV. We should 
also note the quoting of the PDG concerning this reso- 
nance, "There are quite possibly two Poi states in this 
region" . Our results reinforce this hypothesis. 

Finally, in the irirA amplitude for the I = 0, J w = con- 
figuration we find a similar structure at 1740 MeV with 
the full width at half maximum being 20 MeV. We iden- 
tify this peak as the A(1810 — il50/2) resonance, which is 
listed as a three-star l/2 + resonance by the particle data 
group [l|. We would like to draw the attention of the 
reader to the large variation in the peak positions as well 
as the widths reported by different partial wave analyzes 
Q for the A(1810) resonance (the peak position changes 



from 1750 MeV to 1850 MeV and the width from 50-250 
MeV). 

We do not find any states with exotic isospin, I = 2,3. 

We have also investigated the theoretical uncertainties 
of the model. We have already mentioned that P m in and 
the cut off in the three particle loop do not play any rel- 
evant role. In addition, we have checked the sensitivity 
of our results to the change in the two body input pa- 
rameters. We have varied the pion decay constant and 
the two body cut offs, by about 5%, which still guaran- 
tees a fair agreement of our two body cross sections with 
the experimental ones. We find changes in the peak po- 
sitions by less than 5 MeV from each source, or 7 MeV 
when summed in quadrature. This gives us an idea of 
the accuracy of our results. 

The states obtained are not exotic and their quantum 
numbers can be reached with just three quarks. What 
our findings arc telling is that in Nature these three quark 
states unavoidably couple to two mesons and one baryon, 
and, for the states that we have found, the two meson one 
baryon components overcome the weight of the original 
three quarks seed. This particular nature could be tested 
experimentally by means of different reactions, among 
which, the strong three body decay channels and the ra- 
diative decays [23j should play an important role and 
deserve further theoretical and experimental studies. 

We conclude the discussion by emphasizing that all the 
low lying l/2 + £ and A resonances in the PDG [lii], up to 
the 1800 MeV energy region, get dynamically generated 
as two meson-one baryon states in these calculations. In 
addition, we predict the quantum numbers of the £(1560) 
and also find evidence for a l/2 + £ resonance at ~ 1620 
MeV. It is rewarding to see that the widths obtained in 
this work, which correspond to decay into three body 
systems, are smaller than the total ones to which the 
two body decay widths also contribute. There would be 
no contradiction with these two body channels having a 
smaller weight in the resonance wavefunctions, as implic- 
itly assumed in our study, and having a fair contribution 
to the total width, since some of the three body channels 
to which the resonances couple are kincmatically closed 
for decay, and others which are open have a far smaller 
phase space than that available for two body decay chan- 
nels. 
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